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Abstract We present a detailed assessment of the far-infrared confusion noise imposed on measurements with the ISOPHOT 
far-infrared detectors and cameras aboard the ISO satellite. We provide confusion noise values for all measurement configu- 
rations and observing modes of ISOPHOT in the 90 /im < A < 200 fim wavelength range. Based on these results we also give 
estimates for cirrus confusion noise levels at the resolution limits of current and future instruments of infrared space telescopes: 
Spitzer/MIPS, ASTRO-F/FIS and Herschel/PACS. 
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1. Introduction 



Confusion noise is a major limitation in sensitivity and pho- 
tometric accuracy for the measurements performed with the 
far-infrared (FIR) fil ters/detectors of the ISOPHOT instrument 
jLemke e t al.ll 199 61). on-b oard the Infrared S pace Observatory 
(ISO. lKessler et allll996l) . As was shown bv lxiss et alJ ( l200ll 
hereafter Paper I) measurements in the long-wavelength filters 
of the CI 00 camera (90 and 100/^m) were affected roughly 
equally by confusion and instrument noise and measurements 
with the C200 detector were confusion noise limited. In Paper I 
we provided estimates of the sky confusion noise - the sum 
of cirrus confusion noise and fluctuations of the cosmic far- 
infrared background - for a special measurement configuration 
(one target position bracketed by two reference positions, sin- 
gle pixel apertures of the ISOPHOT CI 00 and C200 cameras) 
and for four ISOPHOT filters. However, it is desirable to inves- 
tigate the dependence of the confusion noise on all actual mea- 
surement configurations. This is a fundamental aspect in the 
scientific validation and interpretation of FIR ISOPHOT mea- 
surements. 

Confusion noise predictions for future/current space mis- 
sions working in the far-infrared usually consider the fluctua- 
tions due to the cosmic far-infrared background (CFIRB) only 
(see e.g. Dole et al., 2003, for Spitzer/MIPS, Jeong et al., 2003, 
for ASTRO-F/FIS and Negrello et al., 2004, for Spitzer/MIPS). 
For deep cosmological surveys the cirrus contribution can be 
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minimized by a careful selection of fields with low Galactic 
emission. However, most of the FIR sky is heavily affected by 
this phenomenon. 



The strength of the cirrus confusion noise is believed 
to decrease rapidly with improving spatial resolution (see 
e.g.lGautier et al.Lll992tlMiville-Deschenes et al.L Eool l2003t 



llngalls etaL| ~p004 ) . Ac cording to the formula given by 
,Helou & Beichmanl ( ll99(]l) . the cirrus confusion noise scales as 



(X/D)^-^ with A being the wavelength of the observation and 
D the diameter of the telescope primary mirror In this respect 
the 3.5 m Herschel Space Telescope will be superior to other 
cyrogenic space missions like ISO, ASTRO-F or Spitzer with 
primary mirror diameters D < Im. Although the structure of 
the Galactic ciiTus may change below the ISOPHOT resolution 
limit, our detailed confusion noise study offers the possibility 
to make predictions for other FIR space telescopes, for the first 
time based on observations in the 170 fim range. 



In this paper we present a detailed analysis of the confu- 
sion noise for ISOPHOT measurements performed with the 
P3, ClOO and C200 detectors in various measurement config- 
urations offered by the ISOPHOT Astro nomical Observation 
Templates (see the ISOPHOT Handbook. ILaureiis et al.ll2003l 
for an overview). Based on these results we provide predic- 
tions for the achievable photometric accuracy (cirrus confusion 
noise at the resolution Umit) for the Spitzer/MIPS, ASTRO- 
F/FIS and Herschel/PACS instruments. 
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2. ISOPHOT instrumental set-up and 
observational parameters 

2. 1. Apertures, filters and detector arrays 

Since the strength of the confusion noise is highly wavelength 
dependent (see Paper I), we considered only those ISOPHOT 
filters, where the confusion noise is at least as strong as the 
typical value of the instrument noise. These are the filters with 
central wavelengths A > 90 /^tm. Some filters lack the required 
number of appropriate maps, therefore we restricted our analy- 
sis to the following detector/filter combinations: ClOO: 90 and 
100 urn, C200 : 170a nd 200 fim (see the ISOPHOT Handbook, 
iLaureiis et all l2003l) . Confusion noise analysis has been per- 
formed for single pixels and for the whole detector array field- 
of-view as well. This is necessary, because for some measure- 
ment modes the photometric flux is derived from the summed- 
up fluxes of the array, e.g. for C200 staring, where the source 
is centered on the common corner of the four detector pixels. 
In the case of single pixels the size of the target/reference aper- 
tures are equal to 46" x46"and 92" x92"for the ClOO and C200 
camera, respectively. In the case of the full arrays the size of the 
target/reference apertures are equal to 138"x 138" [3x3 array] 
and 1 84" x 1 84" [2x2 array] , respectively. Although there were 
no suitable "maps" for the P3 detector - 100/^m filter combi- 
nation, this was modeled with the help of maps obtained by the 
ClOO camera in its 100 fim filter The system res ponses of both 
filters are quite similar (see iLaureiis et al.ll2003l Appendix A). 
Five model apertures were constructed using the ClOO detec- 
tor pixel granulation, and corresponding to the 79", 99", 120" 
and 180" circular and to the 127" x 127" rectangular apertures. 
These model "apertures" were generated by 6x6 pixel matri- 
ces on the ClOO maps, with weights for each pixel suitably set 
according to the theoretical footprint value of this pixel relative 
to the centre of the 100 /xm point-spread function (PSF). 

2.2. Covered observing modes 

A detailed description of ISOPHOT's observing modes (AOT 
= Astronomical Obse rving Template) can be found in the 
ISOPHOT Handbook (Laureiis et all l2003l) . Below we sum- 
marize the essentials features for the derivation of the sky noise. 
Table [0 contains a traceability matrix of applicable configura- 
tions per AOT. 

2.2.1. Staring 

Staring observations were mainly performed as on-off mea- 
surements. The on-off distance was freely selectable, but was 
in most cases a few arcminutes. There was the possibility of 
executing a sparse map including several reference positions to 
any of several target positions. Staring observations are com- 
patible with chopping measurements (rectangular or triangular, 
depending on the number of reference positions) in our analy- 
sis. 



2.2.2. Chopping 

Chopping observations could be performed with: 

- one target and one reference position (rectangular chop- 
ping, see Figs. la-c). Chopper throws had to be chosen in 
the range of 90" < 6* < 180"(see TableQ. 

- one target and two reference positions (triangular and saw- 
tooth chopping, see Fig. la-c). Since we do not consider 
any dependence on the direction of the chopper throw rel- 
ative to the target, sawtooth and triangular chopping are 
equivalent in our analysis and are represented by triangular 
chopping. Typical chopper throws are listed in Tabled 

2.2.3. IVIapping 

Confusion noise on maps may be determined in many ways, 
depending on the point-source flux extraction method. This can 
be e.g. similar to some kind of chopping (see above) using sin- 
gle detector pixels as measuring apertures. A typical and ef- 
fective way is the aperture photometry, using a single detector 
pixel or apertures as target aperture and an annulus of pixels 
placed at a specific distance as reference aperture (see Fig. Id). 

2.2.4. IVIini-maps 

Mini-maps represented a special observing mode which was 
basically used for observing point sources with high on-source 
and background redundancy. In mini-map mode the detector 
(ClOO or C200 cameras) moved 'around' the source in a way 
that the source was centered in each detector pixel once during 
the measurement (see Figs, le and f for a schematic representa- 
tion). The confusion noise was calculated according to the steps 
of the detector motion around a (hypothetic) source consider- 
ing the detector pixel with the source centered in as target and 
the other pixels as reference apertures. All pixels were assumed 
to have equal sensitivities. 

2.2.5. Oversampled maps (P32) 

Oversampled maps were performed on a regular grid, com- 
posed of a series of overlapping parallel scans in the space- 
craft y-axis direction. The chopper was used for oversampling 
between individual spacecraft positions along the scan line. 
The same celestial position was observed during several raster 
pointings allowing for elimination of temporal changes in de- 
tector response. Oversamping factors down to 1/3 of the ClOO 
and C200 detector pixel size (15" and 30", respectively) were 
allowed, therefore the internal pixel size of oversampled maps 
is smaller than that of the maps in our database. However, the 
final photometry in an oversampled map can be performed in 
a similar way as in the case of P22 staring raster maps (see 
Sect. 2.2.3). In this case an aperture size corresponding to the 
ones in our sample (Tables 2-5) should be chosen (e.g. ClOO 
pixel or full detector array apertures), i.e. adding up flux values 
of the small map pixels. 
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ISOPHOT AOT 


submode 


detector 




aperture 


number of 
reference positions 


separation 
target-reference 


PHT03 


staring 


P3 




79"-180" 


1 


lf5-10' 




chopping 


P3 




79"-180" 


1-2 


90"-180" 


PHT05 


staring 


P3 




79"-180" 


1 


1^5-10' 


PHT 17/18/19 


staring 


P3 




79"-180" 


>1 


l.'5-l?5 


PHT22 


staring 


ClOO 




138" X 138" 


1 


2:5-10' 




staring 


C200 




184" X 184" 


1 


3'-10' 




chopping 


ClOO 




46"x46" 


1-2 


135"-180" 




chopping 


C200 




184" X 184" 


1 


180" 




mini-map (3 x 3) 


ClOO 




46"x46" 


24 


46"-130" 




mini-map (2 x 2) 


C200 




92"x92" 


8 


92"-130" 


PHT 25 


staring 


ClOO 




138" X 138" 


1 


2^5-10' 




staring 


C200 




184" X 184" 


1 


3'-10' 


PHT 32 


chopping oversampled 


ClOO 


46" 


x46" or 138" X 138" 


>1 


46"-138" 




chopping oversampled 


C200 


92" 


x92" or 184" X 184" 


>1 


92"-184" 


PHT 37/38/39 


staring 


C200 




184" X 184" 


>1 


3'-l?5 



Table 1. Traceability matrix for the configurations of the investigated ClOO, C200 and P3 Astronomical Observing Templates 



3. Confusion noise analysis 

3.1. ISOPHOT maps 

Our database was built from the final maps produced for 
Paper I. Therefore we refer to this paper for a detailed descrip- 
tion of the data reduction. In summary, the data reduction com- 
prised the following main steps: 



- basic data analysis with PIA' V9.0 ( iGabriel et all Il997l) 
from raw data (integration ramps, ERD) to surface bright- 
ness calibrated maps (AAP); 

- flat-fielding using first-quartile normalisation; 

- subtraction of the Zodiacal Light emission; 

- calculation of the instrument noise. 

3.2. Derivation of tine confusion noise 

The confusion noise of a far-infrared ma p is characterized by 
the structure function of Jc*''order (see e.g. lGautier et allll992t 
iHerbstmeier et allll998i for an introduction): 



S{0, k) 



B{x)-\-Y,B{x + e^ 



(1) 



where B is the measured sky brightness, x is the location of 
the target, k is the number of reference apertures, 9_j, is the sep- 
aration vector of the target and the fc*'* reference aperture and 
the average is taken in spatial coordinates over the whole map. 
9_^-s are determined by the actual measurement configuration. 
Measurement configurations investigated in this paper are il- 
lustrated in Fig.n 

The structure noise due to the fluctuations of the sky bright- 
ness and instrument noise, Ngtr, is defined as: 



(2) 



Nstr{e, k) = x/s{0,k) X n 

' PIA is a joint development by ESA Astrophysics Division and the 
ISOPHOT consortium led by the Max-PIanck-Institut fUr Astronomie 
(MPIA), Heidelberg 



where ft is the solid angle of the measuring aperture. As shown 
in Paper I, the relationship between the confusion noise and the 
instrument noise is: 



Ni 



(3) 



where N is the real sky confusion noise and Ninst is the in- 
strument noise characteristic of a specific far-infrared map. We 
calculated using the steps (Eqs. 1-3) described above for 
specific measurement configurations. 

The final confusion noise values were correlated with the 
average surface brightness of the fields, and this relationship 
was fitted by a 3-parameter equation (see Paper I): 



N(g,k,A) 
1 mjy 



= Co(^,k,A) + Ci(e,k,A). 



B 



IMJy sr- 



77(9,k,A) 



(4) 



The Co, Ci and r] parameters are all functions of the num- 
ber (fc) and configuration {&} of reference apertures and the 
wavelength A of the observation. Note that this confusion noise 
is the superposition of extragalactic background and cirrus con- 
fusion noise. 

3.3. Parameter fitting 

The Co, Ci and ?y coefficients were determined by using a rou- 
tine based on standard IDL- fun ctions. The r outine used the 
Levenberg-Marquardt technique jPress et all 1 19921) to solve 
the least-squares problem. 

To perform a completely successful fit for a specific fil- 
ter/configuration combination, it is necessary to have data 
points in the whole surface brightness range. Due to the lack 
of faint fields observed in the ClOO 100 ^m and C200 200 /xm 
filters the Co parameters for these filters cannot be fitted prop- 
erly. Therefore conversions of Cq coefficients have been ap- 
phed from ClOO 90 ^im to 100 and from C200 170 ^im to 
200 /im, as described in Sect. 3.4. 
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Figure 1. Measurement configurations investigated in our analysis. Target and reference apertures are represented by "T" and 
"R" flags, respectively. In our analysis an aperture can be a single detector pixel (P3, ClOO, C200) or a full detector-array (ClOO 
or C200) as well, (a) ClOO chopping configurations. For rectangular chopping only one reference aperture (Rl) and the target 
aperture (T) is used. For triangular chopping both Rl and R2 reference apertures are applied. Both the target and reference 
apertures can be single detector pixels (D) or full detector arrays {D'). (b) The same as (a) but for the C200 detector; (c) The 
same as (a) but for the P3 detector. The P3 detector was mainly used in combination with circular apertures, (d) Circular (annular) 
aperture at 9 distance around the target aperture. The width of the annulus equals the diameter of the target aperture. In practice 
the annular reference aperture is reahzed by a combination of map detector pixels as shown in the figure; (e) Schematic view of 
the ClOO detector mini-map mode (mini-map with 3x3 raster steps). Raster pointings are indicated by the blue dots. The target 
and reference apertures for the first raster pointing are highlighted in colour; (f) Schematic view of the C200 detector mini-map 
mode (mini-map with 2x2 raster steps, see ClOO description) 



Co coefficients for the various P3 apertures were deter- 
mined by simulated confusion noise measurements on syn- 
thetic maps with Poissonian brightness-distribution. In these 
simulated maps the brightness-scaling was arbitrary, and only 
the confusion noise ratio of a single (ClOO) pixel to a simu- 
lated P3 aperture was calculated for a specific configuration. 
Then these ratios were applied to the Co values of the appro- 
priate configurations in Table|3]to obtain the Co coefficients for 
P3 apertures. 



3.4. Coefficients for non-investigated f liters 

For some other filters, which are similary strongly affected by 
confusion noise (ClOO: 105 /im, C200: 120, 150 and 180^m) 
the number of available maps was too small to perform the in- 
vestigation as above. For these filters we use transformed co- 
efficients Co, Ci and of the filters investigated. The transfor- 
mations are done as follows: 

- The Co and Ci coefficients of the confusion noise - sur- 
face brightness relations reflect the spatial structure of 
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Figure!. Demonstration of the accuracy of confu- 
sion noise predictions using transformations of the 
Co, Ci and 77 parameters. The relative uncertainty 
AN/N = \Ntr - N„ieas\/N„ieas IS plotted vcrsus sur- 
face brightness. Ntr is the confusion noise calculated using 
transformed Co, Ci and jy parameters and N„ieas is the confu- 
sion noise calculated using parameters fitted to measurement 
data of this specific wavelength. The displayed cases are the 
90 to 100 /im and 170 to 200/im transformations, respectively. 
The lower and upper limits of the shaded areas reflect the 
variations among different measurement configurations. 

the emission. This is expressed via the spectra l index a 
( iHelou & Beichmarl Il990l: iGautier et all 1 19921) introduc- 
ing a scahng of (Ai/Ao)^"^". 

- Since a « for the extragalactic background, Co scales 
as: Co(Ai,k,0)= Co(Ao,k,0)x(Ai/Ao). 

- In Paper II we derived an average cirrus spectral index of 
{a) ~ -3 for all filters investigated there. We assume here 
that this can be applied to other wavelengths as well. This 
leads to the scaling: Ci(Ai,k,6l) = Ci(Ao,k,6') x(Ai/Ao)^-^. 

- ?7 values show a relatively small scatter, and are indepen- 
dent of the filter or even the detector used. The average val- 
ues for the 90-100 yum and for the 170-200 yum filters (see 
Tables 2.. .6) are 77100 =1.44±0.18 and 77200 =1.58±0.19, re- 
spectively. This implies that either the average values or 
those of individual configurations can be adopted for pre- 
dictions at other wavelengths. Since the ?; values of spe- 
cific configurations may still reflect some individual prop- 
erties, we apply those directly in the transformations, i.e. 
7;(Ai,k,0) = 77(Ao,k,0). 

The usability of these transformed parameters was veri- 
fied for those pairs of filters for which a sufficient number 
of measurements was available. For the ClOO 100 /im filter 
a transformation was done from the ClOO 90 /im filter and 
compared with the 100/im measurements proper, the same 
was done for the C200 170 /im and 200/im filter pair. This 
test was only possible for the Ci and 77 parameters. Both 
the ClOO 100/im and C200 200/im data bases lack low sur- 
face brightness measurements and Co had to be derived from 
the 90 /<m and 170/im value, respectively. The relative un- 



certainty AN/N = \Ntr — N„ieas\/N„ieas, was determined, 
where Nt,. is the confusion noise using transformed parame- 
ters and Nmeas is the confusion noise using parameters fitted 
to measurement data. The results for the 90-to-lOO /tm and the 
170-to-200/im transformations are presented in Fig.|2l The fit- 
ting of the Co, Ci and 77 parameters on the one hand and sys- 
tematic effects like differences in filter bandwidth, changes in 
the steepness or flatness of the sky background spectral energy 
distribution between two filters, etc. on the other hand, intro- 
duce uncertainties. Considering these, our tests proved that the 
conversion scheme can be effectively used to estimate confu- 
sion noise for non-investigated filters. Taking into account the 
dependence of confusion noise on surface brightness and the 
uncertainty of measuring the latter, even an uncertainty of 50% 
provides an acceptable range for the confusion noise estimates. 

4. Results for ISOPHOT filters 

4.1. Tables and data products 



ap. 


con. 


6 


Co 
(mJy) 


Ci 
(mJy) 




P 


R 


92" 


10.5± 2.3 


0.95±0.11 


1.47±0.31 


P 


R 


138" 


8.7± 2.5 


1.97±0.26 


1.27±0.26 


P 


R 


184" 


8.6± 2.7 


2.04±0.26 


1.34±0.26 


P 


R 


230" 


8.3± 2.8 


2.10±0.28 


1.38±0.26 


P 


T 


92" 


9.4± 1.9 


0.63±0.10 


1.47±0.33 


P 


T 


138" 


7.3± 2.0 


2.09±0.30 


1.07±0.24 


P 


T 


184" 


7.9± 2.0 


1.67±0.19 


1.25±0.26 


P 


T 


230" 


6.5± 2.0 


2.43±0.42 


1.17±0.24 


P 


C 


92" 


8.1± 1.8 


0.74±0.14 


1.36±0.30 


P 


C 


138" 


6.3± 1.8 


1.63±0.28 


1.11±0.25 


P 


C 


184" 


6.0± 1.7 


1.56±0.15 


1.23±0.26 


F 


R 


92" 


25.4±14.7 


9.30±2.52 


1.32±0.18 


F 


R 


138" 


28.4±19.3 


13.27±3.79 


1.33±0.17 


F 


R 


184" 


23.8±21.4 


17.35±5.10 


1.32±0.16 


F 


R 


230" 


17.9±21.8 


20.70±6.16 


1.32±0.15 


F 


T 


92" 


28.6±10.2 


2.58±0.43 


1.33±0.26 


F 


T 


138" 


32.6±15.7 


6.46±1.42 


1.31±0.18 


F 


T 


184" 


30.6±16.3 


8.89±3.88 


1.35±0.35 


F 


T 


230" 


25.5±18.4 


12.59±6.11 


1.30±0.34 



Table 2. Fitted coefficients of Eq.Hfor the ClOO camera 90/im 
filter. Abbreviations in the table: "ap." = aperture: P= single 
pixel, F = full array; "con." = configuration for sky reference de- 
terminaion: R = rectangular chopping (one reference position), 
T = triangular chopping (two reference positions), C = annular 
aperture. 



The fitted Cq, Ci and 77 parameters for a specific detector / 
filter / configuration combination are tabulated in Tables|2jQ It 
should be emphasized that the presented confusion noise values 
are 'per beam' and la values. They have to be coiTected with 
the appropriate PSF-fraction of the beam when compared with 
point-source confusion noise values. An example of the func- 
tional behaviour of the confusion noise with surface brightness 
is given in fig. 2 in Paper I. 
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ap. 


con. 


e 


Co 
(mJy) 


Ci 
(mJy) 




P 


R 


92" 


11.6± 2.5 


1.35±0.71 


1.39±0.11 


P 


R 


138" 


9.6± 2.8 


1.47±0.67 


1.44±0.12 


P 


R 


184" 


9.5± 3.0 


1.52±0.72 


1.51±0.13 


P 


R 


230" 


9.2± 3.1 


1.60±0.75 


1.55±0.15 


P 


T 


92" 


10.3± 2.1 


1.14±0.69 


1.29±0.08 


P 


T 


138" 


8.1± 2.2 


1.44±0.70 


1.28±0.12 


P 


T 


184" 


8.7± 2.3 


1.42±0.70 


1.43±0.12 


P 


T 


230" 


7.2± 2.2 


1.57±0.75 


1.58±0.14 


P 


C 


92" 


9.0± 2.0 


1.13±0.60 


1.25±0.08 


P 


C 


138" 


7.0± 2.0 


1.47±0.70 


1.25±0.12 


P 


C 


184" 


6.6± 1.9 


1.47±0.71 


1.34±0.13 


F 


R 


92" 


28.6±16.3 


6.47±2.03 


1.51 ±0.40 


F 


R 


138" 


31.6±21.4 


8.26±2.53 


1.55±0.41 


F 


R 


184" 


25.6±23.7 


9.56±3.34 


1.59±0.42 


F 


R 


230" 


19.9±24.2 


10.33±3.67 


1.61±0.43 


F 


T 


92" 


31.9±11.3 


3.42±0.98 


1.43±0.45 


F 


T 


138" 


36.2±17.5 


3.71±0.57 


1.63±0.70 


F 


T 


184" 


34.0±18.1 


7.38±2.41 


1.47±0.39 


F 


T 


230" 


25.8±20.4 


14.03±4.15 


1.31±0.26 



Table 3. Fitted coefficients of Eq. |3 for the lOO^m filter 
of the CI 00 detector. Co coefficients have been transformed 
from 90 fim resuhs (see Sect. 2.5). For the meaning of labels 
for apertures (ap.) and configurations (con.) see the caption of 
Table 2. 



ap. 


con. 


e 


Co 


Ci 


V 








(mJy) 


(mJy) 




P 


R 


184" 


19.8± 5.4 


0.94±0.42 


1.86±0.21 


P 


R 


276" 


16.6± 6.2 


2.27±0.54 


1.71±0.20 


P 


R 


368" 


12.3±6.8 


3.74±0.58 


1.61±0.19 


P 


R 


460" 


5.6± 7.5 


5.87±0.61 


1.54±0.18 


P 


T 


184" 


10.2± 4.5 


3.01±0.31 


1.47±0.16 


P 


T 


276" 


18.3± 7.0 


1.61±0.22 


1.66±0.12 


P 


T 


368" 


18.9± 6.7 


1.58±0.24 


1.70±0.12 


P 


T 


460" 


14.4± 7.6 


2.06±0.55 


1.66±0.11 


P 


C 


184" 


9.3± 6.7 


3.37±1.01 


1.46±0.17 


P 


C 


276" 


16.4± 8.5 


1.83±0.15 


1.65±0.23 


P 


C 


368" 


13. 8± 6.8 


1.21±0.50 


1.72±0.21 


F 


R 


184" 


49.1 ±2 1.0 


4.72±1.15 


1.78±0.21 


F 


R 


276" 


40.4±24.2 


10.28±1.29 


1.66±0.20 


F 


R 


368" 


20.7±28.2 


17.56±1.38 


1.57±0.19 


F 


R 


460" 


7.5±27.8 


21.15±1.41 


1.59±0.19 


F 


T 


184" 


51.3±17.6 


1.40±0.51 


1.89±0.23 


F 


T 


276" 


59.2±32.4 


2.01±0.54 


1.93±0.24 


F 


T 


368" 


53.6±23.9 


4.55±1.17 


1.76±0.22 


F 


T 


460" 


36.2±25.4 


9.19±1.30 


1.62±0.19 



Table 4. Fitted coefficients of Eq.|3for the 170^m filter of the 
C200 camera. For the meaning of labels for apertures (ap.) and 
configurations (con.) see the caption of Table 2. 



Due to the surface brightness range of the measurements 
used in our analysis the predictions with the coefficients in 
Tables 2-7 are reliable in the I < (B) < 100 MJy sr^^ and in 
the 1 < (B) <200MJysr-i range for the C100/P3 and the 
C200 detectors, respectively. There are no suitable measure- 
ments at high surface brightness and in addition the structure 



ap. 


con. 


9 


Co 
(mJy) 


Ci 
(mJy) 


^1 


P 


R 


184" 


23.3± 6.3 


7.81±2.35 


1.42±0.32 


P 


R 


276" 


19.5±7.3 


10.03±3.23 


1.46±0.32 


P 


R 


368" 


14.5± 8.1 


11.05±3.64 


1.48±0.33 


P 


R 


460" 


6.6± 8.8 


12.34±4.16 


1.49±0.33 


P 


T 


184" 


12.0± 5.2 


7.17±2.10 


1.30±0.29 


P 


T 


276" 


21.6± 8.2 


7.95±2.40 


1.40±0.31 


P 


T 


368" 


22.3± 7.9 


8.22±2.52 


1.45±0.32 


P 


T 


460" 


16.9± 8.9 


9.75±3.14 


1.47±0.33 


P 


C 


184" 


11.0± 7.8 


7.75±2.32 


1.26±0.28 


P 


C 


276" 


19.3±10.0 


8.58±2.65 


1.37±0.30 


P 


C 


368" 


16.3± 8.0 


6.26±1.76 


1.55±0.34 


F 


R 


184" 


57.8±24.7 


18.20±6.62 


1.58±0.35 


F 


R 


276" 


47.6±28.4 


23.71±8.98 


1.61±0.36 


F 


R 


368" 


24.1 ±33.2 


25.68±9.84 


1.64±0.37 


F 


R 


460" 


8.8±32.7 


24.95±9.52 


1.67±0.37 


F 


T 


184" 


60.4±20.7 


14.07±4.88 


1.46±0.32 


F 


T 


276" 


69.9±38.1 


17.36±6.26 


1.55±0.34 


F 


T 


368" 


63.1±28.1 


18.75±6.87 


1.63±0.36 


F 


T 


460" 


42.6±29.8 


14.78±5.36 


1.74±0.41 



Table 5. Fitted coefficients of Eq. I^for the 200 ^m filter of 
the C200 camera. Co coefficients have been transformed from 
the 170 /im results (see Sect. 2.5). For the meaning of labels 
for apertures (ap.) and configurations (con.) see the caption of 
Table 2. 



detector 


filter 


Co 


Ci 


V 






(mJy) 


(mJy) 




ClOO 


90 ^m 


8.1±2.0 


0.74±0.18 


1.36±0.09 


ClOO 


100 ^m 


9.0±3.1 


0.96±0.24 


1.36±0.10 


C200 


170 ^m 


12.2±3.5 


0.65±0.12 


1.63±0.32 


C200 


200 ^m 


14.4±3.8 


2.35±0.98 


1.41±0.19 



Table 6. Fitted coefficients of Eq. HI for the CI 00/90 yum, 
ClOO/lOO^im, C200/170/.tm, C200/200/im detector/filter 
combinations for the simulated mini-map observing mode. Cq 
coefficients for the 100 /^m filter have been transformed from 
the 90 /im results, for the 200 /im filter from the 170 /im re- 
sults. 



of the FIR emission may change significantly above this level 
(see Kiss et al., 2003, hereafter Paper II). Therefore we cannot 
give accurate estimates for very high surface brightness values 
using our current database. 

Based on the results compiled in the tables we con- 
structed all-sky confusion noise maps, one for each detector 
/ filter / configuration combination. Surface brightness val- 
ues for all positions of our 1° resolution grid were derived 
from COBE/DIRBE data with the Zodiacal Light contribu- 
tion removed and interpolated to the ISOPHOT bands, as de- 
scribed in Paper I. The only difference is that the cirrus colour 
temperature was not fixed to 20 K, but was derived from the 
COBE/DIRBE 100, 140 and 240 /im surface brightness values. 
Due to the in general higher uncertainty a lower weight was 
given to the 140 /^m band data. The long wavelength baseline 
of the 100 and 240 /im bands and the COBE/DIRBE surface 
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aper. con. 6 


Co Ci 77 
(mjy) (mjy) 


19" O R 90" 
79" O T 90" 
99" O T 90" 


7.8±1.7 2.25±0.52 1.15±0.19 
6.0±1.2 1.08±0.19 1.13±0.15 
9.0±1.9 0.74±0.12 1.18±0.18 


79" O T 120" 
99" O R 120" 
99" O T 120" 
120"O R 120" 
120"O T 120" 


6.0±1.2 1.93±0.37 1.10±0.18 
10.9±2.3 2.86±0.27 1.16±0.19 

9.0±1.9 1.59±0.22 1.13±0.19 
14.0±3.0 2.46±0.17 1.17±0.18 
11.7±2.4 1.29±0.11 1.15±0.17 


79" O R 180" 
99" O R 180" 

120"O R 180" 
127"n R 180" 

180"O R 180" 


7.8±1.7 3.30±0.37 1.21±0.20 
10.9±2.3 3.37±0.34 1.18±0.19 
14.0±3.0 3.13±0.26 1.17±0.18 
17.1±4.2 3.26±0.51 1.23±0.21 
24.0±5.2 2.33±0.30 1.23±0.20 


120"O R 240" 


14.0±3.0 3.61±0.34 1.18±0.20 



10.0 



Table 7. Fitted coefficients of Eq.|?]for the simulated P3 detec- 
tor - aperture combinations in the 100/.tm fiher. Coefficients 
were derived for only those combinations, which were selected 
for observations . 



brightness accuracies of a few percent provide a final tempera- 
ture accuracy 5T < 0.2 K. 

We calculated the confusion noise from these surface 
brightness values using Eq. 4 and the Cq, Ci and rj param- 
eters corresponding to the actual measurement configuration 
(Tables 2. ..7). The maps are available in the electronic version 
of this paper^ in FITS format. The data format is described in 
Appendix A.l. One advantage of these maps is that they are 
based on a homogeneous all-sky surface brightness calibration. 
In quite a number of cases the determination of the background 
value from the ISOPHOT measurement itself may be less accu- 
rate (e.g. in chopped observations). In this case the confusion 
noise estimated by using Eq. 4 in combination with the mea- 
sured background brightness is unreliable. Our all-sky maps 
provide information on the average confusion noise around the 
target field, averaged over an area of ^1 deg^, which can be 
more reliable, if there are no strong gradients in the background 
brightness over this scale. There are also maps for the filters 
originally not investigated, with the transformations described 
in Sect. 3.4 applied. 

4.2. Noise performance of various observing modes 

In Fig. |3l we compare the confusion noise values obtained 
for various observing modes. We have chosen the rectangular 
chopping with a single pixel aperture and with a separation of 
Omin = 2pixels as reference (denoted 'PR2'). This is compared 
with the results of other configurations (detailed in the figure 
caption). 

The background determination in mini-map mode turned 
out to provide the lowest confusion noise in most cases. The 
second reference position in triangular chopping reduces the 
confusion noise significantly (by 15-50%) compared to rectan- 
gular chopping. Full detector array apertures prove to be 3-5 



cr 

Dl. 




httpV/kisag.konkoly.hu 



PT2 PC2 FR2 FT2 MM 
measurement configuration 

Figure 3. Comparison of the confusion noise values for dif- 
ferent observing modes and measurement configurations. The 
measurement configurations are coded by the labels and the 
separation B is given in detector pixels of the actual detector: 
PT2 - triangular chopping, pixel aperture, 6* = 2; PC2 - circu- 
lar aperture (annulus of pixels), B = 2\ FR2 - rectangular chop- 
ping, detector array aperture, = 2; FT2 - triangular chopping, 
detector array aperture, = 2; MM - mini-map mode. All val- 
ues are given as ratios with regard to the reference configura- 
tion 'PR2': rectangular chopping, single pixel aperture, = 2. 
The ISOPHOT filters are marked by the following symbols: di- 
amond: ClOO 90Atm; triangle: ClOO, 100 /xm; asterisk: C200. 
170 [im\ circle: C200, 200 /im; The vertical bars limited by the 
arrowheads represent the range of ratios for surface brightness 
values of2MJy sr"! < (B) < lOOMJysr"!. 



times more strongly affected by confusion noise than the same 
configuration with single pixels. 

4.3. Noise chiaracterization of individual 
measurements in ttie ISO Data Archive 

The confusion noise can be a severe limit both for the signal-to- 
noise and the photometric accuracy of faint sources. Due to the 
analysis described here it is now possible to determine a robust 
confusion noise estimate for all FIR ISOPHOT measurements 
performed with the AOTs listed in Table 1 and which are com- 
pact source measurements including background reference po- 
sitions. From ISO Data Archive (Kessler et al., 2003) Version 7 
onward the Data Quality Report of these observations flags the 
possible cirrus confusion contamination and an associated cata- 
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logue file gives the confusion noise numbers calculated both for 
the measured ISOPHOT and the COBE/DIRBE background 
values. These can be directly compared with source fluxes ex- 
tracted from the data products. 

5. Cirrus confusion noise predictions for 
infrared space teiescopes 

5. 1. Simulated maps of cirrus structures 

In this analysis we concentrate on the cirrus component of the 
confusion noise and its scaling between instruments with vari- 
ous spatial resolutions. Here we consider neither the contribu- 
tion of the extragalactic background nor the impact by any in- 
strumental effect. The main characteristic of the cirrus emission 
at a specific wavelength is its spatial structure. This is usually 
described by the spectral index, a, of the power spectrum of 
the image, averaged over annuli (see Paper II for a summary). 
With this parameter the power spectrum is P(/) = -Po(///o)"i 
where P{f ) is the power at the spatial frequency / and Pq is 
the power at the reference spatial frequency /q. As was shown 
in Paper II, a may vary with wavelength and surface bright- 
ness. Previous studies (Papers I & II) unraveled the structure 
of the emission for the scale down to the ISOPHOT resolu- 
tion. Extrapolations to better spatial resolutions of future space 
telescopes can be performed using these results and assuming 
that the general behaviour of the spatial structure remains un- 
changed for higher spatial frequencies, i.e. the same fractal di- 
mension / the same spectral index is valid. 

The correct simulation of the entire cirrus structure of a 
specific sky area in the Fourier space would demand the re- 
production of both the power and phase information. However, 
in any autocorrelation analysis like the confusion noise cal- 
culation, only the Fourier power is important, since the two- 
dimensional autocorrelation function is related to the two- 
dimensional power spectrum only: 

c{0)^— p{f')Mf'e)f'df' (5) 

where P{f') is the two-dimensional power spectrum, Jq{x) is 
the circular Bessel-function of the 0*'' kind and 9 is the angu- 
lar separation. The confusion noise N is proportional to C{6) 
or in more complex configurations iV is a linear combination 
of some C(6',)-s. Different FIR instruments sample the power 
spectrum at different spatial frequencies due to their different 
resolving power, resulting in different confusion noise levels 
on the same ( cirrus) structure. Ori ly Pq depends on the surface 
brightness B dOautier et al.[ll992l) . For a power-law type power 
spectrum the ratio of two power levels at two spatial frequen- 
cies is independent of Pq and therefore of B. If one knows the 
Pq-B relation for one specific instrument (e.g. ISOPHOT), it is 
straightforward to derive a P-B relation for another instrument 
using P(/) = Poif / fo)" ■ This is equivalent with a confusion 
noise - surface brightness relation for this other instrument. 

However, for a specific instrument the size of the detec- 
tor pixels and the measurement configuration play an impor- 
tant role as well. Because of this complexity the easiest way to 
compare different instruments is the confusion noise analysis 




Figure 4. Simulated fractal map with a = -3, as seen 
with the spatial resolution of various telescopes/instruments, 
(a) ISOPHOT 170 ^m; (b) ASTRO-F/FIS 170 /im; (c) 
Spitzer/MIPS 160 /zm; (d) Herschel/PACS 175 ^m 



of simulated maps as these would be obtained by these instru- 
ments. This is done in the same way as for the real ISOPHOT 
maps in Sect. 3. 

For ISOPHOT we derived the relation between the cirrus 
confusion noise and the surface brightness of individual fields 
in Sect. 4. From simulated fractal maps the confusion noise 
ratios of ISOPHOT and another instrument can be obtained. In 
this way the cirrus confusion noise of any FIR instrument can 
be connected to the surface brightness. 

To perform the investigation we constructed high- 
resolution (4096x4096 pixels, (X'S pixel size) synthetic images 
for a range of spectral index values (-2.0 >a> -5.0). The gen- 
eration of th e maps is based o n the random recursive fractal 
algoritlim by lElmegreenI lll997ll . 

The high resolution maps were convolved with the beams 
of the actual telescope/detector/filter combinations and then 
sampled according to the size of the detector pixels. The 
Spitzer/MIPS 160/im point spread function is available at the 
Spitzer Science Center website"^. Point spread functions for the 
Herschel/PACS 110 and 175 i-im simulated measuremen ts wer e 
taken from model calculations tOkumura & Longva| l200ll) . 
The ASTRO-F/FIS PSFs were calculated theore tically using 
the latest information on the telescope design jjeong et all 
| 2003i) . The result of this convolution and resampling is shown 
in Fig. El 

The confusion noise analysis was performed in the same 
way as for ISOPHOT observations. The analysis was restricted 
to the triangular chopping configuration, and the separation 6 
was chosen to be equal to the resolution limits of the actual 
telescope/detector/filter combinations. Since confusion noise is 

* http;//ssc. spitzer.caltech.edu/mips/psffits/ 
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especially important for detection of faint point sources, we 
provide la point source confusion noise values after a correc- 
tion for the footprint of the instruments, instead of single-pixel 
("per beam") confusion noise values. 

5.2. Results 

5.2.1 . Conversion coefficients: 

The main parameters of the investigated instruments and the 
conversion coefficients are summarized in Table |8] The con- 
version equation is: 



R 



R^^ X NpHT 



(6) 



where Npur and Nr are the confusion noise values obtained 
by ISOPHOT and another instrument for the same map, re- 
spectively and R^^ is the conversion coefficient. The R^^ fac- 
tor can also be expressed in terms of the Ci and central point 
spread function fraction (fpsf) factors of ISOPHOT and the in- 
vestigated instrument at their resolution limits: 



(7) 



These R^'^ factors of Table 8 can be compared with the sim- 
ple telescope resolution s caling (X/D)^-^ (for a = -3) in the 
iHelou & Beichmanllll990l) formula. For the 85 cm Spitzer tele- 
scope =4.210"^ and for the 3.5m Herschel telescope 
= 1.210^^. An example of a graphical comparison is 
shown in Fig. 5. 

5.2.2. All-sky confusion noise maps: 

Based on the conversion factors derived above we were able to 
produce low spatial resolution all-sky cirrus confusion noise 
prediction maps in a similar way as for the ISOPHOT all- 
sky confusion noise maps. Surface brightness values were de- 
rived from COBE/DIRBE data, following the same scheme as 
in Sect. 4.1. The confusion noise values of the all-sky maps are 
calculated for the configuration PT2 (Tables 2 and 4,6 = 92" for 
ClOO 90 fim and d = 184" for C200 170 /im ) The Co, Ci and 
rj coefficients for this configuration cannot be applied directly, 
since those resulted from the fits to the total surface bright- 
ness, which contains the contribution of the extragalactic back- 
ground as well. Therefore we fitted parameters using a slightly 
modified version of Eq. 4: 



N 



1 mjy 



B(A) 



Bcfirb(A) 



1 MJy sr 



(8) 



BcFiRB is the surface brightness of the co smic far-infrared 
background at a specific wavelength (c.f. |Pei et al.l 1 19991) 
and B — Bcfirb = Bdrr is the cirrus surface brightness. If 
Bcirr =0, then N = Cq, i.e. the confusion noise is purely due 
to cosmic infrared background fluctuations. This indicates that 
C*! • (Bcirr)^ is the pure cirrus confusion noise component. 
We used these CJ' and 77* values, together with the coefficients 
in Table 8 for the creation of the all-sky cirrus confusion noise 



1 0- 



102 




1 



1 00 



(MJysr-1) 



Figures. Demonstration of the effects of the ISOPHOT- 
DIRBE photometric transformation and variable spectral in- 
dex a for the Herschel/PACS 175 i-im filter cirrus confusion 
noise predicions. Solid line: no DIRBE-ISOPHOT transforma- 
tion, constant a; Dashed line: no DIRBE-ISOPHOT transfor- 
mation, variable a; Dash-dotted line: DIRBE-ISOPHOT con- 
version applied, constant a; Dash-triple-dotted line: DIRBE- 
ISOPHOT conversion applied, variable a. The dotted lin e 
shows the prediction according to lHelou & BeichmanI ( Il99(]l) . 

maps. The Cq, C* and if parameters are listed in the headers 
of the FITS files. 

There are two other issues which have to be considered in 
the production of the confusion noise prediction maps: 

1. ) Transformation between the COBE/DIRBE and ISOPHOT 

photometric systems. The all-sky surface brightness maps 
are in the DIRBE photometric system, while the coef- 
ficients are derived in the ISOPHOT system. We used 
the latest available transformation coefficients based on 
the comparison of DIRBE surface brightness values and 
ISOPHOT mini-map background fluxes obtained with 
PIAV10.0/CALG7.0 (Moor, 2003, priv. comm.). The 
transformation equation was: 

B DIRBE = Gam x B^p^r + Offset (9) 

where B^j^gj^ and B^^^ are the COBE/DIRBE and 
ISOPHOT surface brightness values, respectively. The 
transformation coefficients (Gain and Offset) for different 
filters are summarized in Table |9] 

2. ) Constant or variable spectral index (a). As shown in 

Paper II, the spectral index a depends on the surface bright- 
ness of the observed fields and on the observational wave- 
length. The surface brightness dependence of a is approxi- 
mated by: 

B — Bcfirb 



a 



Al X log 



10 



1 MJy sr- 



Ao 



(10) 
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instrument 


filter 




-2.0) 


= - 


-3.0) 




-4.0) 


R^'^ia = - 


-5.0) 


res. limit 


PHT ref . 


ASTRO-F/FIS 


170 ^im 


4.4x10" 


- i 


3.9x10" 


- i 


3.5x10" 


- i 


3.1x10" 


-i 


65'.' 5 


170 /im 


Spitzer/MIPS 


160 nm 


1.6x10" 


-1 


1.2x10" 


-1 


8.7x10" 


-2 


6.5x10" 


-2 


46'.' 5 


170 /im 


Herschel/PACS 


UO/im 


3.0x10" 


-2 


1.3x10" 


-2 


5.3x10" 


-3 


2.4x10" 


-3 


8'.'5 


90 /im 


Herschel/PACS 


175 nm 


1.6x10" 


-2 


7.4x10" 


-3 


3.0x10" 


-3 


1.4x10" 


3 


13'.' 5 


170 )j.m 



Table 8. Conversion coefficients for point sources between ISOPHOT and the specified instruments for cirrus confusion noise 
calculations. The columns of the table contain: 1.) Instrument; 2.) Filter; 3.) Ratio of confusion noise values at the resolution 
limits for a spectral index of a = -2.0; 4.) same for a = -3.0; 5.) same for a = -4.0; 6.) same for a = -5.0; 7.) Resolution limits 
of the cun-ent instrument/filter combination (76"and 143"for the ISOPHOT 90 and 170 /xm filters, respectively); 8.) ISOPHOT 
reference filter 



filter 


Offset 


Gain 




(MJysr^^) 




90Atm 


-1.65±0.06 


l.OOiO.Ol 


100 Atm 


0.09±0.09 


0.85±0.02 


120 pm 


0.62±0.28 


0.95±0.07 


150 nm 


0.28±0.15 


1.02±0.03 


170 fivn 


0.09±0.09 


1.06±0.02 


180 Atm 


-0.50±0.32 


1.15±0.09 


200 ^m 


0.95±0.17 


0.97±0.03 



Table 9. Transformation coefficients between the ISOPHOT 
and the DIRBE photometric systems, according to Eq.|9] 



where B is the average surface brightness of the field af- 
ter the removal of the Zodiacal Light contribution and 
Be FIB.B is the surface brightness of the cosmic far-infrared 
background. The coefficients are: Ag =-1.67 ±0.47 and 
Al =-1.57±0.38 for long wavelength filters (170/200 ^tm) 
of the C200 detector (see Paper II). For shorter wavelength 
filters (CI 00, 90 & 100 fim) this relation cannot be properly 
derived, therefore we do not use variable spectral indices 
for A < llQiim. 

Due to these considerations two/four maps can be produced for 
each instrument: 

i) No DIRBE - ISOPHOT transformation, constant a = -3 

ii) No DIRBE - ISOPHOT transformation, 
variable a (only for A > 120yLtm) 

iii) DIRBE - ISOPHOT transformation applied, 
constant a = -3 

iv) DIRBE - ISOPHOT transformation apphed, 
variable a (only for A > 120 fim) 

The detailed technical description of these all-sky fits files is 
given in Appendix A. 2. 

The effect of these corrections is demonstrated in Fig.|5]for 
the Herschel/PACS 175 ^lm filter. The effect of the DIRBE- 
ISOPHOT transformation applied or not is not very significant. 
In the contrary, the variable a modifies the slope of the Ndrr- 
Bcirr relation, leading to a factor of ^^5 lower cirrus confusion 
noise at high surface brightness values. 

An example comparing the capabilities of four in- 
struments (ISOPHOT nO^im, ASTRO-F/FIS 170 /im, 
Spitzer/MIPS 160/im, Herschel/PACS 175 /im) is given in 
Fig.El 



6. Summary 

In this paper we present the results of a detailed investiga- 
tion of the dependence of the sky confusion noise and in par- 
ticular the cirrus confusion noise on measurement configura- 
tions for the long-wavelength (A > 90 /im) observations with 
ISOPHOT. Tables 2-7 together with Eq. 4 provide an easy 
tool to estimate confusion noise values. Based on these re- 
sults and transformations for non-investigated filters we con- 
structed all-sky confusion noise maps (stored as FITS files at: 
^ttp://kisag. konkoly.hu/ISO/toolbox.html) for all possible mea- 
surement configurations of the P3 100 /im, CI 00 90, 100 and 
105 ^im and C200 120, 150, 170, 180 and 200 ^im filters. These 
files can efficiently be used to estimate the confusion noise at 
any specific sky position even if the background surface bright- 
ness cannot be properly estimated from the measurement itself. 
From Version 7 onward the ISO Data Archive provides con- 
fusion noise values for individual ISOPHOT measurements, 
based on our numbers. 

The results of the ISOPHOT confusion noise analysis and 
the utilization of simulated fractal maps allowed us to calcu- 
late cirrus confusion noise value ratios of ISOPHOT and other 
far-infrared spaceborn instruments. Using these values, all-sky 
maps with cirrus confusion noise estimates for point sources 
were constructed. This was done for the resolution limits of 
detectors of Spitzer/MIPS, ASTRO-F/FIS and Herschel/PACS. 
These maps can be used for the preparation of FIR observations 
with future space telescopes indicating sensitivity limits due 
to cirrus confusion noise. However, for passively cooled large 
telescopes (like HERSCHEL) other no ise components l ike the 
thermal telescope background (see e. g. lOkumural bOoH) or the 
extragalactic background will play a more dominant role. 
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